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Abstract 

We study the fluctuation properties of transition intensities applying a recently proposed gen- 
eralization of the random matrix theory, which is based on Beck and Cohen's superstatistics. We 
obtain an analytic expression for the distribution of the reduced transition probabilities that ap- 
plies to systems undergoing a transition out of chaos. The obtained distribution fits the results of 
a previous nuclear shell model calculations for some electromagnetic transitions that deviate from 
the Porter-Thomas distribution. It agrees with the experimental reduced transition probabilities 
for the 26 Al nucleus better than the commonly used x 2 distribution. 
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I. INTRODUCTION 



Random matrix theory (RMT) is believed to describe quantal systems whose classical 

nnn 

counterpart has a chaotic dynamics |2J, |3| . RMT the m atrix elements of the Hamiltonian 
in some basis are replaced with random numbers. The theory is based on two main assump- 
tions: (i) the matrix elements are independent identically-distributed random variables, and 
(ii) their distribution is invariant under unitary transformations. These lead to a Gaussian 
probability density distribution for the matrix elements, P (H) oc exp [— r/Ti (H^HM . With 
these assumptions, RMT presents a satisfactory description for numerous chaotic systems. 
Agreement with RMT is now considered to be a signature of chaos in the quantum system. 
For time-reversal-invariant systems, the appropriate form of random matrix theory is the 
Gaussian orthogonal ensemble (GOE) ; that is the form which will mainly be considered in 
this paper. 

For most systems, however, the phase space is partitioned into regular and chaotic do- 
mains. These systems are known as mixed systems. Attempts to generalize RMT to describe 
such mixed systems are numerous; for a review please see Most of these attempts are 
based on constructing ensembles of random matrices whose elements are independent but 
not identically distributed. Thus, the resulting expressions are not invariant under base 
transformation. The first work in this direction is due to Rosenzweig and Porter 4|. They 
model the Hamiltonian of the mixed system by a superposition of a diagonal matrix of 
random elements having the same variance and a matrix drawn from a GOE. Therefore, 
the variances of the diagonal elements of the total Hamiltonian are different from those of 
the off-diagonal ones, unlike the GOE Hamiltonian in which the variances of diagonal el- 
ements are twice those of the off-diagonal ones. Hussein and Pato |^ used the maximum 
entropy principle to construct "deformed" random- matrix ensembles by imposing different 
constraints for the diagonal and off-diagonal elements. This approach has been successfully 
applied to the case of metal-insulator transition 61. A recent review of the deformed en- 

Q 

semble is given in [2J. Ensembles of band random matrices, whose entries are equal to zero 
outside a band of limited width along the principal diagonal, have often been used to model 
mixed systems 0, j], E| . 

The past decade has witnessed a considerable interest devoted to the possible generaliza- 
tion of statistical mechanics. Much work in this direction followed Tsallis seminal paper 1(]J . 
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Tsallis introduced a non-extensive entropy, which depends on a positive parameter q known 
as the entropic index. The standard Shannon entropy is recovered for q — 1. Applications 
of the Tsallis formalism covered a wide class of phenomena; for a review please see, e.g. [ill ]. 

ude systems with mixed regular-chaotic 



)pliea to include systen 

Q H Q, H Q Q 



However, the constraints of 



Recently, the formalism has been ap plie d to inc 
dynamics in the framework of RMT 
normalization and existence of an expectation value for Tr (H^H) set up an upper limit for 
the entropic index q beyond which the involved integrals diverge. This restricts the validity 
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of the non-extensive RMT to a limited range near the chaotic phase 

Another extension of statistical mechanics is provided by the formalism of superstatistics 
(statistics of a statistics), recently proposed by Beck and Cohen Superstatistics arises as 
weighted averages of ordinary statistics (the Boltzmann factor) due to fluctuations of one or 
more intensive parameter (e.g. the inverse temperature). It includes Tsallis' non-extensive 
statistics, for q > 1, as a special case in which the inverse temperature has a ^-distributions. 
This formalism has been applied to model a mixed system within the framework of RMT 



in Ref. 
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20]. 



The joint matrix element distribution was represented as an average over 



exp [— r?Tr (if^if)] with respect to the parameter r]. The different choices of parameter 
distribution, which had been studied in Beck and Cohen's paper were considered in 



The parameter distribution has also been estimated 



maximum entropy, as done by Sattin 



by applying the principle of 
2l\ . Explicit analytical results were obtained for the 



level density and the nearest neighbor-spacing distributions. 

Matrix elements of transition operator probe the system's wave functions so that their 
statistical fluctuations provide additional information. In chaotic systems, the reduced tran- 
sition probabilities follow the Porter-Thomas distribution j^]. This is a ^-distribution of 
one degree of freedom. As the system becomes more regular, the transition probabilities 
deviate from the Porte-Thomas distribution. To account for these deviations, Alhassid and 
Nov 0S e lsky S n gg e st ed t na« the transition widths in odxed system may b e anai yZ ed in 
terms of a ^-distribution of a lower degree of freedom. The latter distribution does not fit 
well the empirical distributions but consists with the observed number of weak transitions 
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). The distributions 
28l | and 30 P Q expressed as 



as compared with the Porter-Thomas distribution (see, e.g. 
of experimental reduced transition probabilities B in 26 A1 
functions of the log-B have peaks at logi? < while all the y 2 distributions are peaked at 
\ogB = 0. We show in the present paper that the superstatistical RMT provides us with 
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a more suitable generalization of the Porter-Thomas distribution. In Section II we briefly 
review the concept of superstatistics and the necessary generalization required to express 
the characteristics of the spectrum of a mixed system into an ensemble of chaotic spectra 
with different local mean level density. The evolution of the reduced transition-intensity 
distribution during the stochastic transition induced by increasing the local-density fluctua- 
tions is considered in Section III. Section IV demonstrates the quality of fit achieved by the 
obtained transition-intensity distribution by comparing its prediction with the results of a 
shell-model calculation by Hamoudi et al. The conclusion of this work is formulated 

in Section 5. 



II. SUPERSTATISTICAL RMT 



To start with, we briefly review the superstatistics concept as introduced by Beck and 
Cohen Q|. Consider a non-equilibrium system with spatiotemporal fluctuations of the 
inverse temperature (3. Locally, i.e. in spatial regions (cells) where j3 is approximately 
constant, the system may be described by a canonical ensemble in which the distribution 
function is given by the Boltzmann factor e~@ E , where E is an effective energy in each 
cell. In the long-term run, the system is described by an average over the fluctuating f3. 
The system is thus characterized by a convolution of two statistics, and hence the name 
"superstatistics". One statistics is given by the Boltzmann factor and the other one by the 
probability distribution f{(3) of j3 in the various cells. One obtains Tsallis' statistics when f3 
has a x 2 distribution, but this is not the only possible choice. Beck and Cohe n g ive several 
possible examples of functions which are possible candidates for /(/?)• Sattin [2l| suggested 
that, lacking any further information, the most probable realization of f(f3) will be the one 
that maximizes the Shannon entropy. Namely this version of superstatistics formalism will 
now be applied to RMT. 



A. Joint distribution of the matrix-elements 



Gaussian random-matrix ensembles have several common features with the canonical 
ensembles. In RMT, the square of a matrix element plays the role of energy of a molecule 
in a gas. When the matrix elements are statistically identical, one expects them to become 
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distributed as the Boltzmann's. One obtains a Gaussian probability density distribution of 
the matrix elements 

P G (H) oc exp [-r/Tr (H^H)] (1) 

by extremizing the Shannon entropy [lj subjected to the constraints of normalization and 
existence of the expectation value of Tr(wH) . Here Tr means trace and is the Hermitian 
cojugate of H. The quantity Tr (H^H) plays the role of the effective energy of the system, 
while the role of the inverse temperature (3 is played by i], being twice the inverse of the 
matrix-element variance. 

Our main assumption is that Beck and Cohen's superstatistics provides a suitable de- 
scription for systems with mixed regular-chaotic dynamics. We consider the spectrum of 
a mixed system as made up of many smaller cells that are temporarily in a chaotic phase. 
Each cell is large enough to obey the statistical requirements of RMT but has a different 
distribution parameter rj associated with it, according to a probability density f(rj). Con- 
sequently, the superstatistical random-matrix ensemble that describes the mixed system is 
a mixture of Gaussian ensembles. Its matrix-element joint probability density distributions 
obtained by integrating distributions of the form in Eq. (1) over all positive values of rj with 
a statistical weight f(r]), 

where Z(rj) = J exp [— rjTr (H^H)] drj. Here we use the "B type superstatistics" The 
distribution in Eq. (2) is isotropic in the matrix-element space. Relations analogous to Eq. 
(1) can also be written for the joint distribution of eigenvalues as well as any other statistic 
that is obtained from it by integration over some of the eigenvalues, such as the nearest- 
neighbor-spacing distribution and the level number variance. The distribution f{j]) has to 
be normalizable, to have at least a finite first moment and to be reduces a delta function as 
the system becomes fully chaotic. 

An analogous ensemble made of a superposition of random-matrix ensembles has recently 



been considered by Muttalib and Klauber 



271 ]. These authors have been seeking for gen- 



eralizations of Gaussian random- matrix ensembles, with the probability distributions P(H) 
that are functions of the single variable Tr (H^H') like the distribution (2) that follows here 
from the concept of superstatistics. However, they work not directly with the distributions 
P(H) themselves, but with the associated characteristic functions defined as the Fourier 



transforms 

C(T) = J exp [iTr (T*H)] P(H)dH. (3) 

They prove, among other things, that if C(T) is a function of Tr (T^T) only, then the 
most general C(T), valid for random-matrix ensembles of arbitrarily large dimension can be 
represented as 

C(T) = J f(b) exp [-bTr (TtT)] db. (4) 

The inverse Fourier transformation of C(T) then leads to an expression for P{H) similar to 
the one in Eq. (2). We consider their result as a justification of using Eq. (2) for ensembles 
of matrices of dimensions N — > oo. 



B. Marginal distribution for a single matrix-element 

Unlike the Gaussian random-matrix ensembles, the superstatistical ensemble has cor- 
related matrix elements. This can clearly be seen by the fact that the joint distribution 
function defined by Eq. (2) does not factorize into a product of distributions of the individ- 
ual matrix elements. However, it is not difficult to obtain a marginal distribution for each 
of the individual matrix elements. 

We shall confine our consideration for the GOE; the generalization to the other symmetry 
universalities is straightforward. In this case, 

Tr (H*H) = Tr(H 2 ) = £ H\ k + 2 £ H%, (5) 

k k<l 

Integrating the joint distribution P(H) over all the matrix elements except one, say Hif, we 
obtain 

roc 

p if (H if )= /Mv^exp {-2r]H? f } drj. (6) 
Jo 

The parameter r] essentially defines the energy scale of the individual ensembles, whose 
superposition compose the superstatistical ensemble. We therefore assume that the distri- 
bution (6) will hold for the superstatistical distribution of any physical quantity having the 
dimension of energy, which is represented as a Gaussian random variable in the case of GOE. 

We note that the distribution function of superstatistical ensemble depends on the matrix 
elements through Tr(H 2 ) which is base invariant. In other words, the fuction P(H) is 
invariant under rotation in the spcae of matrix elements. Therefore the distributions Pif 
have the same form for all off-diagonal matrix elements of H. 
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C. Parameter distribution 

Following Sattin (2^, we use the principle of maximum entropy to evaluate f{f])- Lacking 
a detailed information about the mechanism causing the deviation from the prediction of 
RMT, the most probable realization of f{rf) will be the one that extremizes the Shannon 
entropy 

S = - f( V )\nf( V )d V (7) 
Jo 

with the following constraints: 

Constraint 1. The major parameter of RMT is 77 defined in Eq. (1). Superstatistics 
was introduced in Eq. (2) by allowing r\ to fluctuate around a fixed mean value (rj). This 
implies the existence of this mean value 



(V) = / f(v)vdri- (8) 







Constraint 2. The fluctuation properties are usually defined for unfolded spectra, which 
have a unit mean level spacing. The mean level density is proportional to the inverse square 
root of rj. We thus require the existence of the integral 



f( V )V 1/2 dv = 1- (9) 
Therefore, the most probable f(r]) extremizes the functional 

/*oo poo /»oo 

F = - / f( v ) hif(rj)drj - X, / f( V )rjd V - A 2 / f( V )v~ 1/2 dv, (10) 
Jo Jo Jo 



where Ai and A 2 are Lagrange multipliers. As a result, we obtain 

1/2 



f{r]) = Cexp 



(11) 



^0 V V 

where a and 770 are parameters, which can be expressed in terms of the Lagrange multipliers 
Ai and A 2 , and C is a normalization constant. The latter is given by 

(a*| 0,1,1)- 1 ' 

Here (x\ bi, 62, ^2) is a Meijer's G-function 0, [y]; see also the appendix of Ref. j^J. 
The parameter distribution f(rj) is peaked at rjQ and tends to a delta function as a — > 00. 
The value of t] will be fixed in the next section while the parameter a will be considered as 
the tuning parameter for the stochastic transition. 



A parameter distribution analogous to f{rf) is obtained in [20j, where the variable rj 
is replaced by the local mean level density. This distribution is used in (^(J to derive 
expressions for the level density distribution, the nearest-neighbor spacing distribution, and 
the two-level correlation function for spectra of superstatistical ensembles. The expression 
obtained for the level-density distribution has a finite value at the center of the spectrum 
for arbitrarily large ensemble dimension N, and thus satisfies a necessary condition on f(rj), 
which is required by Muttalib and Klauber [27 1. 

III. TRANSITION-INTENSITY DISTRIBUTION 

The probability Bif of a transition from the initial configuration \i) to the final configu- 
ration \f) is given by 

B lf = \W lf \\ (13) 

where 

W if = (f\0\i) (14) 

is the square of the transition operator O in a special basis. In a chaotic system, the 
eigenstates \i) and |/) are believed to be very complicated. If the operator O conserves 
time reversibility, the matrix elements Wif are real. For a chaotic system, it is reasonable 
to assume that Wif are identically-distributed Gaussian random variable. This entails that 
the transition intensities can be represented by a random variable that takes the values 

mf = t§A- (15) 



where (By) is a suitably defined local average value [24j . and has a Porter-Thomas 
distribution 

PMv) = J^er™. (16) 
V ny 

The parameter r\ is defined by the requirement that (y) = 1, and is equal to 1/2. A more 



elaborate derivation of the Porter-Thomas distribution is given by Barbosa et al. [25 1. 

We now derive the superstatistical generalization of the Porter- Thomas distribution. For 
this purpose, we assume that the matrix elements Wif are distributed according to Eq. (6). 
The parameter rj in Eq. (16) is no more considered as a constant but allowed to fluctuate 



S 



according to the distribution f(rj). The superstatistical transition intensity distribution is 
then given by 

r°° ~ nf _ 

-Psuperstati S tical(|/) = / f(v)\ e VV d V ( 17 ) 

Jo V 

Substituting Eq. (11) for f(r)) and integrating over i], we obtain 

P a ^(a 2 (a + w)|0,il) riR , 

Superstatistical [V) ~ ^ ^ + a/r]o )3/2 G 30 ^ Q> 1 ^ ^ l A «J 

The parameter ?7q is determined from the requirement that (y) = 1, which yields 

a gg (o 3 |0,0, |) 
"° 2G»(a'|0,i,l)' 1 ' 

Replacing Meijer's G-function by its large- argument asymptotic expression 

Gl'l (Z \h, b 2 , b 3 ) ~ ^ Z (bi+^+b 3 -l)/3 exp (_ 3 ^l/3) (20) 

v3 

one can easily show that -Psu P orstatisticai(l/) is reduced to the Porter-Thomas distribution as 
the parameter a — > oo. 

Several independent results with different models, have suggested that transition 
strengths in a chaotic system follow a x 2 distribution 

PAv>») = ^kjg) yl,/2 ' le ~ vv ' (21) 

with v — 1 (porter-Thomas) degrees of freedom, the transition strengths in a less chaotic 
system ay 2 distribution with a number of degrees of freedom less than one. Alhassid and 
Levine |32| introduced this distribution using maximum-entropy arguments. Several studies 



of electromagnetic transition intensities in nuclei have also been performed 
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24. 
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29J; each of these has suggested that the x 2 distribution with v < 1 is appropriate for 
relatively regular systems. 

Experimental data for transition intensities range over several orders of magnitude. It is 
often more convenient to consider a logarithmic variable as the argument. The probability 
density function in terms of log a y is 

F(\og a y)=y\naF(y). (22) 

We compare in Fig. 1 the evolution of -Psu P crstatisticai(logio2/) during the transition from 
chaotic to regular dynamics by varying the tuning parameter a from oo (GOE) to 10 -6 



(almost regular) with a corresponding evolution of P x 2 (logioy, v) where v varies between 
1 and 0.1. Of special interest is the fact that the maximum of P x 2(logioy, v) occurs at 
logioy = for any value of v. This property does not hold for Psuperstatisticai(logio2/)- The 
peak of superstatistical distribution for less chaotic systems occurs at logioy < and moves 
towards lower values as the parameter a decreases. We show in the next section that this is 
indeed the behavior of physical systems. 



IV. DATA ANALYSIS 



The purpose of this section is to show that the proposed superstatistical distribution 
succeeds in the situations where the x 2 distribution fails. We show this by using results 
from two works, which examine the effect of a transition from chaos to integrability on 
gamma-ray reduced transition probabilities. 

The first work is done by Hamoudi, Nazmitdinov and Alhassid j^. They calculated 
the electric quadrupole (E2) and magnetic dipole (Ml) transition intensities among the 
isospin T = 0, 1 states of nuclei with mass number 60. They applied the interacting shell 
model with realistic interaction for pf shell nuclei with a 56 Ni core. It is found that the 
B(E2) transitions are well described by a GOE (Porter-Thomas distribution). However, the 
statistics for the B(M1) transitions is sensitive to T z . The M2 transition operator consists 
of an isoscalar and isovector components. The T z — 1 nuclei, in which both components 
contribute, exhibit a Porter-Thomas distribution. In the meanwhile, a significant deviation 
from the GOE statistics for the T z = nuclei, where the matrix elements are purely isoscalar 
and relatively weak j33j]. 

We analyzed the reduced Ml transition intensities for both the T z — 1 60 Co nuclei and 
T z = 60 Zn calculated by Hamoudi et al. 2(| nuclei using the superstatistical transition 
intensity distribution in Eq. (14). These authors sampled a large number of matrix elements 
for each transition operator, which is equal to 56 2 = 3136 and 66 2 = 4356. Figure 2 compares 
the results of calculations using Eq. (14) with the numerical results of Hamoudi et al. [3| 
as well as the "best-fit" \ 2 distribution deduced by these authors. The figure clearly shows 
the advantage of the superstatistical distribution proposed here over the x 2 distribution, at 
least for this numerical experiment. 

The second work that we consider here is that of Adams, Mitchell and Shriner (3]. They 
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collected approximately 1500 experimental reduced electromagnetic transition strengths be- 
tween the excited states of the nucleus 26 Al. Their data involve levels with isospin T = and 
T = 1 between the ground state and the excitation energy of 8.067 MeV. Figure 3 compares 
these experimental data with results of calculations using Eq. (14) with a = 1.24 as well as 
the "best-fit" x 2 distribution with a parameter v slightly greater than 1. The figure again 
shows the advantage of the superstatistical distribution over the x 2 distribution, although 
the agreement with the data is not as good as in the cases shown in Fig. 2. The experi- 
mental histogram is mostly higher than the theoretical curves especially in the peak region, 
although the data was normalized to to 0.83 in order to approximately take care of the 
upper and lower detection thresholds. The percentage of undetected transitions may have 
been underestimated because its estimation was based on the Porter-Thomas distribution. 



V. CONCLUSION 



The eigenstates of a chaotic system are extended and cover the whole domain of classically 
permitted motion randomly, but uniformly. They overlap substantially, as manifested by 
level repulsion. There are no preferred eigenstate; the states are statistically equivalent. As 
a result, the matrix elements of transition operators in any basis are independent and have a 
Gaussian distribution, which leads to the Porter-Thomas distribution for reduced transition 
intensities. Coming out of the chaotic phase, the extended eigenstates become less and 
less homogeneous in space. Different eigenstates become localized in different places and 
the matrix elements that couple different pairs are no more statistically equal. The matrix 
elements will no more have the same variance; one has to allow each of them to have its 
own variance. But this will dramatically increase the number of parameters of the theory. 
The proposed superstatistical approach solves this problem by treating all of the matrix 
elements as having a common variance, not fixed but fluctuating. One then expresses the 
probability density of transition intensities as an average of Porter-Thomas distributions 
with different mean intensities. The principle of maximum entropy is used to estimate the 
inverse-mean-intensity distribution. The resulting transition-intensity distribution is found 
to agree with the results of shell model calculation as well as with experimental data better 
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that the x 2 distribution, which is often used for this purpose. 
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Figure Caption 

Figure 1. (Color on line) Evolution of the superstatistical distribution 
-Psu P crstatisticai(logio2/) and the x 2 distirbutions P x 2(log w y,v) during the transition from 
chaotic to regular dynamics. The solid curves, labeled as PT, refer to the Porter-Thomas 
distribution. 

Fig. 2. (Color on line) Nuclear shell-model Ml transition intensities in A = 60, cal- 

n 

culated by Hamoudi et al. |26j, (histograms) compared with the superstatistical distribu- 
tion (14) with parameters a = 1.499, 0.064 and 0.030, respectively (solid curves) and the 
^distribution (17) with parameters v — 1, 0.64 and 0.34, respectively (dashed curves). 
Fig. 3. (Color on line) The distribution of experimental reduced transition probabilities 



in 26A1 from Ref. 



24j (histogram) compared with the superstatistical distribution (14) with 



parameter a = 1.24 (solid curve) and the ^distribution (17) with parameters v = 1.04 
(dashed curve). 
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